One-dimensional orbital fluctuations and the exotic magnetic properties of YVO3 
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Starting from the Mott insulator picture for cubic vanadates, we derive and investigate the model 
of superexchange interactions between V 3+ ions, with nearly degenerate ti g orbitals occupied by two 
electrons each. The superexchange interactions are strongly frustrated and demonstrate a strong 
interrelation between possible types of magnetic and orbital order. We elucidate the prominent role 
played by fluctuations of yz and x z orbitals which generate ferromagnetic superexchange interactions 
even in the absence of Hund's exchange. In this limit we find orbital valence bond state which 
is replaced either by C-type antiferromagnetic order with weak G-type orbital order at increasing 
Hund's exchange, or instead by G-type antiferromagnetic order when the lattice distortions stabilize 
G-type orbital order. Both phases are observed in YVO3 and we argue that a dimerized G-type 
antiferromagnetic phase with stronger and weaker FM bonds alternating along the c axis may be 
stabilized by large spin-orbital entropy at finite temperature. This suggests a scenario which explains 
the origin of the exotic C-AF order observed in YVO3 in the regime of intermediate temperatures 
and allows one to specify the necessary ingredients of a more complete future theory. 
[Published in: Phys. Rev. B 75, 184434 (2007).] 
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ORBITAL DEGREES OF FREEDOM 



Transition metal oxides with the perovskite struc- 
ture display a large variety of properties such as high- 
temperature superconductivity and colossal magnetore- 
sistance. Their magnetic properties are also quite di- 
verse, with antiferromagnetic (AF), disordered, or ferro- 
magnetic (FM) phases in different doping regimes, be- 
ing the subject of particularly active research in the last 
decade. 1 ' 2 Although certain universal principles can be 
formulated, these complex magnetic properties depend 
on the actual filling of 3d orbitals of transition metal ions, 
and have to be studied in detail for each family of com- 
pounds separately. Rich and complex behavior in doped 
systems is found as moving charges can dress by spin or 
orbital excitations.- The undoped compounds are some- 
what simpler as their properties are dominated by large 
on-site Coulomb interactions oc U, responsible for their 
Mott-Hubbard (or charge-transfer) insulating behavior, 
with the effective low-energy magnetic interactions of su- 
perexchange type. While such interactions are AF and 
nonfrustrated on a cubic lattice for nondegenerate or- 
bitals, they have a very nontrivial structure when degen- 
erate 3c? orbitals are partly occupied, as pointed out by 
Kugel and Khomskii on the example of e g systems long 
ago 4 In such cases the orbital degrees of freedom have 
to be considered on equal footing with electron spins^ 
which leads to the so-called spin-orbital superexchange 
modelsj£i£&£ describing the low-energy physics and the 
partial sum rules in the optical spectroscopy*^ 

An intriguing feature of the spin-orbital models is the 
strong frustration of the superexchange interactions on 



a cubic (perovskite) lattice which was recognized as the 
origin of enhanced quantum effects in transition metal 
oxidesiii For purely electronic models this frustration 
might even lead to the collapse of long range order in 
particular parameter regimes, but usually this does not 
happen and the fluctuations are partly suppressed ei- 
ther by the order-out-of-disorder mechanism^ or by the 
coupling to the lattice distortions induced by the Jahn- 
Teller (JT) effect. In the cuprates and manganites both 
the superexchange and the JT interactions support each 
other ; 13 : 1 and such systems undergo usually structural 
transitions. Nevertheless, even below the structural tran- 
sition the spin and orbital degrees of freedom are cou- 
pled, leading to characteristic changes of the orbital order 
(00) at magnetic transitions and to new composite spin- 
orbital excitations, when both spin and orbital excitation 
occurs simultaneously^ 

The importance of the orbital degrees of freedom was 
realized in the theory of magnetism already in the sev- 
enties. Next to e g systems^ model Hamiltonians with 
twofold degeneracy and diagonal hopping, 17 and the real- 
istic effecrive Hamiltonian for ti g electrons in V2O3 were 
studied^ Actually, the superexchange interactions for 
partly filled t2 g orbitals are different and even more fasci- 
nating than those for e g systems. As realized first for the 
d 1 configuration in cubic titanates^ the quantum effects 
are here even stronger than in the e g systems (cuprates 
or manganites), as the JT coupling is weak and the or- 
bitals may form the coherent orbital liquid state observed 
in a Mott insulator LaTiO3, 20 As a result of this quan- 
tum behavior and common spin-orbital fluctuations, the 
classical Goodenough-Kanamori rules are violated in ti g 
systems in some cases ^ 
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The quantum effects are equally important in vana- 
dium compounds with V 3+ ions in the d 2 configuration, 
realized in V2O3 and in cubic compounds: LaVC>3 and 
YVO3. The metal- insulator transition in V2O3 is stud- 
ied for quite a long timeji but more realistic superex- 
change models were introduced only after the experimen- 
tal evidence of the 00 which occurs below the magnetic 
transition^ The first spin-orbital model for V2O3 was 
assuming a picture of molecular bonds which saturated 
one t 2g electron per V 3+ site and thus used s = 1/2 
spins However, one decade ago it was realized that 
Hund's exchange Jh is large^ and the superexchangc 
interactions couple instead 5 = 1 spins of different V 3+ 
ions. A complete superexchange model with spin and or- 
bital degrees of freedom in V2O3 was derived only a few 
years ago by Di Matteo, Perkins and Natolii2i 

Also in cubic vanadates evidence increases that the or- 
bital degrees of freedom couple to the magnetic order and 
play an important and highly nontrivial role. In LaV03 
the C-type AF (C-AF) phase [with FM chains along c 
axis which stagger within ab planes] is stable below Neel 
temperature TV — 143 K, followed by a weak structural 
transition at T s ~ 141 K^^^^^^ Remarkably, the 
magnetic order parameter in the C-AF phase of LaV03 
is strongly reduced to ~ 1.3/isi^ As the spin quantum 
fluctuations are smaller than in the G-AF phase and are 
unlikely to decrease the order parameter by more than 
6% for 5=1 spins^i the observed large reduction of 
the magnetic moments suggests that some other quan- 
tum effects which originate from the orbital degeneracy 
dominate in this phase of cubic vanadates. 

The situation is very different and even more puzzling 
in YV0^22&&21& — this compound has G-type AF 
order (staggered in all three directions, called below G- 
AF) at low temperatures T < TV2, while the magnetic 
order changes in the first order magnetic transition at 
Tn2 = 77 K to G-AF structure which remains stable 
up to Tjvi — 116 K. The magnetic transition at TV2 is 
particularly surprising as the staggered moments change 
their direction from approximately parallel to the c axis 
in the G-AF phase to lying almost within the ab planes in 
the G-AF phase, with some small alternating G-AF-like 
component. 3 - 3 " In addition, the magnetization is strongly 
reduced at T > Xat 2 , being only close to 1.0/zb in the 
G-AF phased and the magnetic exchange constants J a b 
and I J c \ are there much lower than those found in the low- 
temperature G-AF phased Even more surprising is the 
observed gap in the spin wave spectrum which suggests 
an exotic dimerized structure with alternating stronger 
and weaker FM exchange constants along c axisi 36 i 37 In 
addition, recent Raman experiments 3 ^ suggest that the 
short-range orbital fluctuations of the orbital G-type oc- 
cur in this intermediate G-AF phase in addition to the 
alternating orbital (AO) G-type (G-AO) order, and make 
it thus quite different from the one observed in LaV03. 
We also note that the competition between G-AF and 
G-AF phase is a common feature of a few vanadate com- 
pounds with low atomic radiii 3 ^ 



The electronic structure calculations gave valuable in- 
formation about the possible charge distribution over the 
t2 g orbitals in YV03i 38 i 39 Large on-site Coulomb inter- 
action U prevents double occuppancy of d orbitals — it 
is implemented in the calculations using the local den- 
sity approximation (LDA) within the so-called LDA+/7 
method^ The commonly accepted picture is that the 
xy orbitals are occupied by one electron, while the sec- 
ond one occupies either yz or xz orbital. The lattice 
distortions in YVO3 are larger in the low-temperature 
phase and suggest G-AO order. Above TV2 the dis- 
tortions decrease and are compatible with a weak G- 
type AO (G-AO) order i^l Theoretical analysis within the 
charge-transfer model has shown that both phases are in- 
deed energetically closed and one may thus expect that 
small changes of the thermodynamic potential around 
T/V2 could induce a first order phase transition. 

In this paper we study the magnetic properties of cu- 
bic vanadates with a spin-orbital model derived for vana- 
dates some time ago42 This model applies to Mott in- 
sulators with transition metal ions with partly filled t 2g 
orbitals in either d 2 or d 4 configuration. Therefore, this 
model was recently used to analyze the magnetic struc- 
ture of monolayer ruthenates^ 3 - In the context of vana- 
dates we have already shown before that the orbital fluc- 
tuations play a prominent role in this model and amplify 
the FM coupling along the c axis, providing a microscopic 
explanation of the observed G-AF order in LaV03. In 
fact, FM interactions induced by Hund's exchange cx Jh 
alone are typically much weaker than the AF ones, and 
would not be sufficient to explain why the FM interac- 
tions are even stronger than the AF ones in the high 
temperature G-AF phase of YVO3. 

Here we will concentrate on the exotic magnetic prop- 
erties of YVO3 and address several open questions moti- 
vated by the observed magnetic properties, in particular 
why: (i) the spin exchange interactions are so different 
in G-type and G-type AF phases of YVO3, (ii) the mag- 
netic transition at TV 2 takes place, (iii) the order param- 
eter (S z ) in the G-AF is so strongly reduced, and finally, 
(iv) the dimerization along the c axis, observed in the 
G-AF phase in the intermediate regime of temperature 
T/V2 < T < T/vi, takes place. A careful discussion of 
these questions in the context of the microscopic model 
will lead us to a scenario for the exotic magnetic proper- 
ties of the intermediate temperature phase of YVO3 con- 
sistently explained within a dimerized G-AF order stable 
only at finite temperature, and characterized by reduced 
exchange interactions. At the same time, we will argue 
that further theoretical studies are necessary in order to 
explain all the observed properties. 

The paper is organized as follows. In Sec. |TT]we present 
the spin-orbital model for cubic vanadates. It is derived 
from the degenerate Hubbard model (Sec. Ill A[) and con- 
tains superexchange interactions supplemented by orbital 
interactions induced by the lattice (Sec. IIIB|) . Next we 
introduce the possible types of classical order in Sec. IIII1 
emphasizing first the tendency towards one-dimensional 
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(ID) orbital fluctuations (Sec. IIII A[) , and next compar- 
ing their classical energies (Sec. IIIIBI) . The effective 
exchange interactions in different magnetic phases are 
evaluated in Sec. IIV Al For the magnetic phases stable 
in different regimes of parameters we derive spin (Sec. 
IIV Bp and orbital (Sec. IIV C|l excitations, which serve 
next to calculate the quantum corrections to the energy 
and lead to the phase diagram of the model at T = of 
Sec. IPTDl 

Using the above background information we propose 
a scenario for the magnetic phase transition at T/v2 in 
YVO3 in Sec. [V] The unique instability of the ID spin- 
orbital chain (Sec. IV A[) comes here together with the 
reduction of the magnetic exchange constants by orbital 
fluctuations (Sec. IV Bp to stabilize the dimerized C-AF 
phase at temperature T > T^2, as we show by analyzing 
the spin and orbital entropy contributions to the free en- 
ergy (Sec. IV Cp . In Sec. IVII wc summarize the results and 
present general conclusions. The paper includes two ap- 
pendices which present the derivation of the spin-orbital 
model for cubic vanadates (Appendix [A]) , and the calcu- 
lation of spin and orbital excitations, as well as the av- 
erage order parameters, and intersite (spin and orbital) 
correlations at finite temperature in the dimerized C-AF 
phase (Appendix iBjl . 

II. SPIN-ORBITAL MODEL FOR CUBIC 
VANADATES 

A. Degenerate Hubbard model for t<z a electrons 

We consider a realistic degenerate Hubbard model for 
3d electrons of V 3+ ions in cubic vanadates, with partly 
filled t2 g orbitals that are energetically favored over e g 
orbitals by the octahedral field. Thereby, we neglect 
small lattice distortions and the tilting of VC>6 octahedra. 
Therefore, the e g orbitals do not couple to t 2g orbitals by 
the hopping processes and play no role in the magnetic 
properties we address below. In such an (idealized) per- 
ovskite structure V 3+ ions occupy the cubic lattice, and 
the hopping elements between active t^g orbitals are the 
same in all three cubic directions. The model Hamilto- 
nian, 

H = H t + H ci + H int , (2.1) 

includes the kinetic energy H t , the orbital splittings in- 
duced by the crystal field H ci , and the on-site electron- 
electron interactions H lri t- The kinetic energy is de- 
scribed by the effective hopping element t between two 
V 3+ ions which originates from two hopping processes via 
the 2pn oxygen orbital along each Mn-O-Mn bond. Its 
value can in principle be derived from the charge-transfer 
model; 3 -^ 3 . and one expects t — i? pd j 'A ~ 0.2 eV. A more 
accurate estimation from the theory is not possible at the 
moment, so we will have to rely on experimental informa- 
tion from neutrobn scattering concerning the magnetic 
exchange constants in YVO3. 



The kinetic energy is given by: 

flt = -*E E ( c C<w + 4/-<v) . ( 2 - 2 ) 

where d\ va are electron creation operators, and the sum- 
mation runs over the bonds {ij) \\j along three cubic axes, 
7 = a, b, c. As observed before ; 19 ! 42 only two out of three 
t2 g orbitals, labelled by ^(7), are active along each bond 
{ij) and contribute to the kinetic energy (|2.2p . while the 
third orbital lies in the plane perpendicular to the 7 axis 
and the hopping via the 2p^ oxygen is forbidden by sym- 
metry. This motivates a convenient notation used below, 

\a) = \yz), \b) = \xz), \c) = \xy), (2.3) 

with the inactive orbital along a given cubic direction 7, 
labelled by its index as I7). 

The electron-electron interactions are described by the 
on-site terms ^ 

- 2J H E S ifl -S iu + J H E d irt d inl d i»l d ivp 

(2.4) 

with U and Jh standing for the intraorbital Coulomb and 
on-site Hund's exchange interaction, respectively, using 
the notation of Kanamori4£ Each pair of orbitals {/x, 1/} 
is included only once in the respective interaction terms 
with summations over [i < v. The Hamiltonian (|2.4p 
describes rigorously the multiplet structure of d? and 
d 3 ions within the t^g subspace 4 ^ and is rotationally in- 
variant in the orbital spaced 4 - More precisely, the on-site 
Coulomb interactions depend on three Racah parameters 
{A, B, C}, and for t2 g orbitals one finds, 

U = A + AB + 3C*, J H = 35 + C, (2.5) 

The Coulomb and exchange element, U and Jh, can be 
thus obtained using the spectroscopic information about 
the Racah parameters for V 2+ ions in the excited states: 
A = 3.54 eV, B = 0.095 eV, and C = 0.354 eV, as given 
by Zaanen and Sawatzky4i With these parameters one 
finds U = 5.0 eV and J H = 0.64 eV. 

The Coulomb element U is therefore sufficiently large 
compared to t ~ 0.2 eV (i.e., U ^> t) to use the second 
order perturbation theory in which the charge fluctua- 
tions dfd'j ^ d 3 dj are suppressed, and the d electrons are 
localized in tf g configurations of a Mott insulator (The 
interaction parameters for V 3+ ions have similar values 
to those of V 2+ ones) . We use this picture as a starting 
point for our analysis and assume that two electrons are 
localized at each U 3+ ion i, satisfying a local constraint 
(at site i) for the total electron density, 

n t = n la + n ib + n ic = 2, (2.6) 
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FIG. 1: An artist view of the energy splittings between t2 g 
orbitals in YVO3 in different temperature regimes. The or- 
bital splitting A which occurs below the structural transition 
at T s and persists in the C-AF phase favors the occupied xy 
orbitals, but allows also for weak orbital fluctuations. Such 
fluctuations are quenched in the G-AF phase at T < Tn2- 



where rii a = ^^riiaa, etcetera. Two electrons at ev- 
ery site are in the high-spin 3 T 2 triplet (S = 1) state, 
stabilized by Hund's exchange Jh- As t <C Jh, the ki- 
netic energy H t can only contribute in virtual processes 
which are responsible for the superexchange interactions 
derived below in Sec. Ill Bl 

The third term in Eq. (|2.1[) stands for the orbital en- 
ergies in crystal field induced by the structural transition 
at T s ~ which lifts the de generacy of three ti g 

orbitals and breaks the cubic symmetry in the orbital 
space. We write the crystal field term H c { as follows, 



S = 1 spin forms at each ion, and one finds a large de- 
generacy 9 N of the ground state, where N is the number 
of sites, as every spin component (S z = 1,0,-1) is al- 
lowed, and a hole may occupy either orbital: |a), \b) or 
|c). This large degeneracy is, however, removed by the 
effective interactions between each pair of nearest neigh- 
bor ions {«, j}, which originate from virtual transitions to 
the excited states due to charge dfd 2 ^ g? 3 dj excitations, 
generated in each case by a single hopping of a ti g elec- 
tron. In the realistic regime of parameters such processes 
may be treated perturbatively, and one arrives in second 
order perturbation theory at an effective superexchange 
Hamiltonian of Ref. UJ| — the details of the derivation 
are explained in Appendix [A"l 

The superexchange interactions between two S = 1 
spins at sites i and j arise from virtual excitations 
d 2 d 2 — > tifdj along the concerned bond (ij), promoted 
by the hopping t which couples pairs of identical active 
t2 g orbitals. A single hopping process generates a <i 3 con- 
figuration, cither with three different orbitals occupied by 
a single electron each, or with a double occupancy in one 
of the two active orbitals (see Fig. Therefore, the df 
excited state may be either a high-spin 4 ^2 state, or one 
of three low-spin states: 2 E, 2 T\ or 2 T 2 with energies^ 8 - 
U — 3Jh, U and U + 2J#, as shown in Fig. 1 of Ref. d. 
This perturbative consideration leads to the spin-orbital 
superexchange model for cubic vanadates, 



Hj 



J 



Si ■ s, + 1 J, 



?(7) 



K 



(7) 



(2.9) 



(2.7) 



with the energy scale given by the superexchange con- 
stant, 



4f 
U ' 



(2.10) 



with electron energies for orbital fi at site i. In agree- 
ment with the results of band structure calculations , 38 ! 39 
and with an idealized but suggested by the local distor- 
tions and thus commonly accepted picture^ we assume The spin interactions oc Si ■ Sj obey the SU(2) symme- 

(7) 



that the xy orbitals are favored below the structural tran- 
sition, while the remaining yz and xz orbitals are nearly 
degenerate, i.e., e c < e a , and Sb — e a , leading to 



1, 



1. 



(2.8) 



i.e., c orbitals are 'condensed' and the other two represent 
the remaining t2 g orbital degree of freedom at every site 
(see Fig. [TJ). Although in principle the orbital energies 
Ei^ could change at the magnetic transition at Tjv2 and 
further stabilize G-AF phase at low temperature, we will 
ignore small corrections which would result from this ef- 
fect in the derivation of the superexchange, and consider 
only generic features of the spin-orbital model that could 
be responsible for the experimental situation. 



B. Superexchange model for vanadates 

Consider first the atomic limit, i.e., the system of V 3+ 
ions in d 2 configuration at t = 0. In the ground state 



try. In contrast, the orbital interaction operators J n 

and K^p involve only two active t2 g orbitals on each in- 
dividual bond (ij) || 7 (7 = a, b, c) which contribute to 
the virtual excitations, so they have a lower (cubic) sym- 
metry. These operators take the form: 



J 



(7) 



(1 + 2777-1) \ T i -T j + jn i n j 



7?r 3 ( fi ® fj + -n i n j ) - -7/Ti(7ij + rij) 
1 \ ( ^ 1 



1(7) 



(2.11) 



(7) 



1 



-(1 +Wi){n, ,+rij) 



(2.12) 



and have a rich structure which originates from the pro- 
jections of the df excited states on the respective eigen- 
states of V 2+ ion, as explained in Appendix [A"l 
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FIG. 2: Virtual charge excitations dfdj — * d^d^ — t djd 2 
within a bond (ij) along c axis, which contribute to the su- 
perexchange in cubic vanadates. Orbital fluctuations which 
support the FM superexchange occur when different active 
orbitals a and b are occupied at both sites, as in cases (a) 
and (b). If the same orbitals are occupied at both sites, e.g. 
the orbital a as in case (c), the superexchange is AF — then 
a double occupancy of the occupied (and active) orbital is 
generated in the excited state, which next dissociates to a 
configuration with either: (c) the same orbital occupancies, 
or (d) with interchanged occupied orbitals at sites i and j. 



First of all, the interactions J^J and K^p depend on 
Hund's exchange splittings in the multiplet structure of 
a V 2+ ion in local d 3 configuration (shown in Fig. 1 of 
Ref. 0) via the exchange parameter, 



V = 



J_h_ 

U 



(2.13) 



and the respective coefficients r\ and r$ in Eqs. (|2. 1 1|) 
and (|2.12p are (for convenience we use here the same 
notation as in Ref. 0): 



1 



n = 



1-37?' 



•'•3 



1 + 2?/ 



(2.14) 



They correspond to the excitation spectrum in dfd 2 f± 
d?dj charge transitions (Fig. In the present case of 
cubic vanadates one finds 4 ^. j) ~ 0.13, which we will take 
as a representative value for YVO3. 

The pseudospin (orbital) operators = t~ , rf } 
for pseudospin r = 1/2 in Eqs. (|2.1ip and (|2.12p are 
defined in the subspace spanned by two orbital flavors 



which are active along a given direction 7. For instance, 
the virtual transitions which generate the superexchange 
interactions follow from the electron hopping between the 
pairs of active a and b orbitals along the bond (ij) \\ c 
axis (see Fig. [2]), and these operators are defined by Eqs. 
(|A6j) . while the number of active electrons at site i is 

(c) 

n \ = n ia + n ib- ^ * s important to realize that although 
the pseudospin flavor is conserved in each individual hop- 
ping processes, the off-diagonal elements of the Coulomb 
interaction oc generate transitions between the com- 
ponents of the excited states, as shown in Fig. [^d). 49 
Therefore, next to the usual scalar products, 



1 



» 3 ' a i 3 



- I f; • f; - -ii; 11 : j = {n ia n ja + aj&fbjoj-) + (a <-» b), 

(2-15) 

we also find in the orbital operators and 'orbital 
fluctuating' terms 



(c) 



= ( n ta n ja + 4 b i a \ b i) + (« ^ 6 )' 



(2.16) 

where (a <-> b) stands for the terms with interchanged 
a and b orbitals. Unlike in the Heiscnberg model, the 
interactions oc t^t^ = ajb^bj in Eq. (|2.16|) induce 
similar orbital flips at both sites. Such terms have the 
form 



(n 



i) — 



1 



't?t1 



(2.17) 



and lead to the nonconservation of the total pseudospin 
quantum number and are thus responsible for further en- 
hancement of orbital quantum fluctuations on the bonds 
with both orbitals active (in this case along c axis). In 
contrast, the bonds in ab planes are classical as there 
analogous terms cannot contribute when the c orbitals 
have condensed. This demonstrates that the breaking 
of symmetry in the orbital space, such as given by Eqs. 
(|2.8|) , will have severe consequences for magnetism. 

The complete microscopic model we consider in the 
following Sections, 



Tt = Ttj + Ti., 



orb ; 



(2.18) 



includes as well effective orbital interactions induced by 
the oxygen distortions. When the VC>6 octahedra dis- 
tort at a second magnetic transition at T, 



N2- 



32,35 



intersite 



interactions which help to order yz and xz orbitals, oc- 
cupied by one electron et every site, are induced. They 
are of two types — the GdFe03-type distortions favor 
repeated orbitals along the c axis, while the AO order 
in ab planes is favored by weak JT effect. Therefore, in 
addition to the superexchange (|2.9[) we introduce two ef- 
fective orbital interactions {V C ,V Q } as the last term of 
the effective Hamiltonian (|2.18p . 



7~Lorb — 



(i3)\\c 



V a 



2^ T i T 3 



(2.19) 



<jj)||afc 
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where the orbital pseudospin operator t[ at site i is de- 
fined by Eq. (j A7|) . With the present sign convention 
both parameters are positive (V c > and V a > 0) and 
induce the C-type AO (C-AO) order, as observed in the 
G-AF phase at T < TV 2 . For convenience we express 
the orbital interactions in Hoib (|2.19p in the units of the 
superexchange constant J, and introduce dimensionless 
parameters: 



(2.20) 



V a V c 

v a = - v c = -, 

which describe the model given by Eq. (|2.18[) , in addition 
to Hund's exchange parameter r\. 



III. 



TYPES OF MAGNETIC ORDER 



Orbital singlets at Jh 







In order to understand the possible symmetry breaking 
in the cubic vanadates, consider first the superexchange 
interactions in the Jh —> limit: 



Si ■ Sj + 1 ) In ■ f. 



1 

-nm, 
4 1 j 



(t) 



(3.1) 

where a constant energy of —2 J per V 3+ ion is neglected. 
It is straighforward to understand why the interactions 
at Jh turn out to have the same structure as in 
LaTi03fi£ where for spins s = 1/2 of Ti 3+ ions one finds 
instead the spin part 4(s*j • Sj + k). In fact, in the limit 
of Jh — > the superexchange interactions follow entirely 
from the Pauli principle, as the multiplet structure of 
excited states collapses to a single degenerate level and 
the spin interactions cx Si ■ Sj due to the high-spin 4 A-i 
and low-spin 2 E states, which involve d 3 {abc} configura- 
tions, cancel each other (see Appendix lAf. This suggests 
that the superexchange interactions might all be AF in 
the limit of Jh — ► 0, as in e g systemsiiiii In fact, in e g 
systems only one directional orbital is active along the 
bond, two electrons occupying these orbitals form an in- 
traorbital spin singlet, which maximizes the energy gain 
for the AF superexchange. 

However, there is an important difference between the 
e g (with one hole per site) and ti g (with one or two 
electrons per site) systems, which may be best realized 
by considering a single bond (ij) in one cubic direction. 
Two active t2 g orbitals along this bond open a new pos- 
sibility — if both orbitals are singly occupied, an or- 
bital singlet gives here FM superexchange, even in the 
absence of Hund's exchange Jh^ For the present filling 
of n = 2 electrons per site, and if n c = 1, such a reso- 
nance on a bond is possible only along one out of three 
cubic directions^ — the orbital singlets and uncorrelated 
bonds alternate along the c axis [Fig. [3ja)]. In analogy 
with spin systems^ this state can be called an orbital 
valence bond (OVB) stated This possibility was also in- 
dependently pointed out by Shen, Xie and Zhang^ 3 - who 
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FIG. 3: (Color online) Schematic picture of the classical 
phases with magnetic and orbital order for ni C — 1 in ac 
plane: (a) OVB phase, with alternating strong FM bonds 
stabilized by orbital singlets represented by double lines, and 
weak AF bonds (with disordered a/b orbitals); (b) C-AF spin 
order accompanied by G-AO order; and (c) G-AF spin order 
accompanied by G-AO order with repeated either a or 6 or- 
bitals along c axis. In cases (b) and (c) spins and orbitals 
alternate along b direction (not shown). These latter states 
follow the Goodenough-Kanamori rulea^ 1 ^ and are analyzed 
below for YVO s . 



obtained the OVB state as the most stable solution of 
the present Hamiltonian (12. 9ft in the regime of small 77 
and for large S limit. 

The OVB state implies an unconventional type of mag- 
netic order. At 77 = the exchange constants along the 
c axis are given by 



Jc{v = 0) 



1 



(3.2) 



When the orbital singlets form and contribute to the en- 
ergy with (fi ■ Tj)( c ' — —j, they maximize the FM ex- 
change on these bonds [see Sec. IIV A| . and stabilize there 
effective 5 = 2 spin states. Between them one finds dis- 



ordered orbitals, i.e. 



(n ■ Tj)( c ' = 0, so the magnetic 
exchange interactions on these bonds are much weaker 
and are in fact AF due to the static term (n,*n,-)^ = 1 



in Eq. (|3.2[) . The interactions within the ab planes are 
also AF but somewhat stronger — they follow from the 
conventional (Pauli principle) mechanism which operates 
as well in the absence of orbital degeneracy, with intraor- 
bital singlets generated by the nearest neighbor hopping 
between sites with singly occupied c orbitals (12.811 . As- 



suming disordered orbitals one finds (fi • fj)W> = and 
{ninj)( ah } = I for the bonds in ab planes. The resulting 
magnetic order which coexists with the orbital singlets is 
shown in Fig. (Ha). 



B. Magnetic and orbital order at finite Jh 

Let us analyze the possible types of coexisting mag- 
netic and orbital order of the full effective Hamiltonian 
given by Eq. (|2.18p for finite Hund's exchange Jh, which 
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includes the effective orbital interactions {14,14} with 
the lattice. Motivated by the experimental situation in 
YVO3, we assume that the c orbitals have condensed, 
so the constraints given by Eqs. (|2.8[) arc fulfilled. For 
this case we consider possible classical phases and their 
energies. A more complete analysis which includes the 
quantum corrections due to spin and orbital excitations is 
presented in Sec. IIVDI here we discuss only a qualitative 
picture when spin quantum fluctuations are neglected. 

At r\ = the lowest energy is obtained when the orbital 
fluctuations are fully developed at every second bond 
along the c axis in the OVB state ; 36 ' 37 ' 53 as discussed 
in Sec. MI Al The classical energy of this phase per site 
is obtained assuming the classical values for intersite spin 
correlations: (Si ■ Sj) = ±1 on the FM/AF bonds. It in- 
cludes the orbital fluctuation energy gained on the orbital 
singlet bonds and is given by: 



p(o) _ _ T 
^OVB — J 



11. . 1 

-n + ~r){9ri - llr 3 ) - -v c 



(3.3) 



Here and below we neglect a constant nonmagnetic term 
— 2J. Except for the orbital singlets, the orbital interac- 
tions in 7i or b do not contribute as the {a, 6} orbitals are 
disordered on all other bonds [Fig. Ela)]. 

An alternative AF state, realized at larger values of 
r)£L is obtained when the (negative) orbital correlations 
along the c axis are uniform, and all the bonds exhibit 
FM exchange. As the spin interactions remain AF in ab 
planes, these interactions lead to the C-AF phase shown 
schematically in Fig. [Hb). A straightforward estimate 
of the classical energy of this phase, 



4° 3 = J 



l\(c) 



ri{fi ■ Tj + -J + r](ri + r 3 )(n ia n ja ) 



(>') 



i](2ri - r 3 ) - v c + 2v a 



(3.4) 



is obtained again taking the classical spin correlations: 
(g t ■ Sj)^ = 1 and (St ■ Sj)^ = -1. It depends on 
the orbital correlations (f ■ fj). Taking fully disordered 
ID orbital chain with (f • fj) = —0.4431, as for the AF 
Heisenberg spin chain^ one finds a crossover from the 
OVB to the C-AF phase at ry — 0.064. We improve this 
naive estimate of the transition in Sec. IIV D[ where we 
evaluate the quantum corrections due to spin excitations 
in both phases. 

Unlike in e g systems ; 13 ' 1 the orbital interactions in- 
duced by the lattice (]2.19[) compete with the superex- 
change (|2.9[) in the present effective spin-orbital model 
(|2.18[) and stabilize the G-AF phase at sufficiently large 
orbital interaction V c £^ The classical energy of this 
phase, 



E, 



(0) 
G 



J 



^ 1 

V(ri - r 3 ) + -v c + 



1 



(3.5) 



is lowered by the energy — | J(v c + 2v a ) gained per site 
when the C-AO order shown in Fig. (3Jc) sets in. In fact, 



the C-AO order enforces here the G-AF phase, show- 
ing a close interrelation of spin and orbital intersite cor- 
relations, known in the literature as the Goodenough- 
Kanamori rules 



IV. SPIN AND ORBITAL EXCITATIONS 

A. Effective exchange interations 

In order to analyze the spin and orbital excitations, 
we follow the usual approach in mean field theory^ and 
decouple spin and orbital operators in Eq. (|2.9[) . Note 
that this approach is satisfactory below the spin order- 
ing temperature Tjvi, as then the spin fluctuations are 
quenched and the spin and orbital degrees of freedom 
may be disentangled^ while for T > Tjvi composite 
spin-orbital excitations need to be considered. This pro- 
cedure leads to the effective spin exchange constants J c 
and J a b, as given in Refs. [9lfl0l. 



Jc 



Jab — 



If. . 

^2 J \n r i - i r i - vn - v r 3) 

j J 1 - 777-1 - i]r 3 



(c) 



+(n - wi - ^3)(rfr/ + - 



(4.1) 



(4.2) 



They depend on the orbital correlations, (fi ■ fj) and 
( T i T j) a l° n g c axis, and (t?t?) in ab planes, which have 
to be determined from the full superexchange model 
given by Eq. (|2. 18[) . i.e., in presence of orbital inter- 
actions promoted by the lattice. Below we specify the 
effective exchange interactions for three possible phases 
shown in Fig. [3] 

At low 77 one expects that the OVB state with alter- 
nating FM and AF bonds along c axis is stable [Fig. 
|3fa)]. On the bonds occupied by orbital singlets, with 



\W = 
\-« ■ j 1 

FM exchange 



and (TfTf) (c) = 



Jcl^—^Jnii + v), 



one finds strong 



(4.3) 



which is further enhanced with increasing 77 and soon 
becomes the dominating magnetic interaction, see Fig. 
|4|a). In contrast, for the bonds connecting singlets the 



resulting AF exchange interactions 



1 



J c i = z J{l-il(2r 1 +r 3 )} 1 



(4.4) 



decrease with increasing 77. These exchange interactions 
are much weaker than the AF ones in the ab planes, 



•Jnh 



1G 



J[l-r){n + r 3 )}, 



(4.5) 
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FIG. 4: (Color online) Exchange interactions J c and J a (, as 
functions of Hund's exchange r\ as obtained for: (a) OVB 
phase with alternating strong FM J c \ (solid line) and weak 
AF J C 2 (dashed-dotted line) exchange interaction; (b) C-AF 
phase with (weak) G-type OO; and (c) G-AF phase stabilized 
by the orbital interactions {14,14} which induce the G-AO 
order. FM (AF) exchange interactions along c axis in G-AF 
(G-AF) phase are shown by solid lines, while AF interactions 
in ab planes are shown by dashed lines. 



constants in the C-AF phase, and J^ b , can be found 
from the orbital excitations in the ID disordered orbital 
chain, and we provide analytic expressions to evaluate 
them in Sec. IIV CI They allow one to determine the 
(weak) OO parameter (r z ) and the intersite orbital cor- 
relations which appear in Eqs. (|4.ip and (|4.2|) . Here we 
present the result of the numerical calculation, see Fig. 
S[b). The FM exchange J~ is finite already at r\ = 
due to the a/b orbital fluctuations^ and is further en- 
hanced by increasing splitting between the high-spin and 
low-spin excitations when Hund's exchange 77 increases. 
At the same time, the AF exchange interaction in ab 
planes decreases. 

Finally, we consider G-AF phase realized in YVO3 at 
low temperature T < T/v2- A classical state with robust 
C-AO order has been proposed for this phased as shown 
in Fig. [3fc) . We have verified that the quantum correc- 
tions to the OO parameter (r z ) are indeed negligible by 
considering the orbital waves for such a classical C-AO 
phase, see Sec. IIVCI so one finds indeed rather simple 
expressions for the AF exchange constants along the c 
axis and in ab planes: 



J ab 



4 

~J[l-r)( ri +r 3 )}. 



(4.7) 
(4.8) 



Both above coupling constants decrease with increasing 
Hund's exchange, and the anisotropy between J c and J a b 
is gradually enhanced (Fig. 2}. 



Spin wave excitations 



in the entire allowed regime of 77, as the latter interactions 
are supported by the excitations of doubly occupied con- 
figurations in c orbitals. One finds that the OVB state 
with alternating FM and AF bonds is destroyed at a crit- 
ical value of 77, 



1 



2r x + ra 



0.188, 



(4.6) 



where the weaker AF bond J C 2 collapses and changes 
its sign, see Fig. |4}a). In reality, it turns out that the 
orbital singlets are destabilized even much faster as a 
better energy is obtained when the spins reorient to FM 
order and the C-AF phase with uniform disordered (or 
weakly ordered) {a, 6} orbitals along the c axis takes over, 
as we show below. 

The simplest possible approach to the C-AF phase is 
to assume that the coupling to the lattice dominates and 
stabilizes the G-AO order, as shown in Fig. ^b). Such 
a robust OO would lead to classical values of orbital cor- 
relations in G-type phases, with (fi • fj + = 0, and 
( T i T j) = 0- However, it was recently pointed out 9 
that this situation does not occur in LaV03, and instead 
one has to consider fluctuating orbitals. The exchange 



The spin waves in different phases can be derived us- 
ing the linear spin wave (LSW) theor y 55 i 56 In the present 
case of S = 1 spins this approach gives reliable results 
also for the OVB phase, in constrast to the linear orbital 
wave (LOW) theory for r = 1/2 pseudospins which can- 
not be applied to the OVB phase as the a and b orbitals 
are there disordered. For the AF phases with two (or 
four) sublattices and the classical AF order (Sf) = ±S 
considered here, we first rotate the spin operators on the 
sites occupied by down spins (with (Sf) — —S) by angle 
7r with respect to spin x axis, which leads to the canonical 
transformation: 



S 



-Sf 



-Sf 



(4.9) 



Next we write the equations of motion for the 
spin operators and apply standard Holstein-Primakoff 
transformation^ from spin operators to boson operators 
(here S = 1), 



S 



2Sa, 



s: 



2Sa 



Sf = S-a]a t . (4.10) 



The respective boson problem is easily diagonalized by 
employing first the Fourier transformation and next a 
Bogoliubov transformation in the momentum space k. 



9 




FIG. 5: (Color online) Spin-wave dispersion ojk (full lines) 
and orbital triplet excitation energy f2k (dashed line), as ob- 
tained for the OVB phase along high symmetry directions in 
the Brillouin zone at r\ = and V c = 0. For the present 
parameters one finds the following values of spin exchange 
constants: J ab = 0.3125J, J c i = -0.250J and J c2 = 0.125J. 
The high symmetry points are: V = (0,0,0), M — (n,n,0), 
R = (7r,7r,7r), Z = (0,0, tt). 



Following the above procedure, one finds the spin- wave 
dispersion in the OVB phase, 



(k) = { [A.J ab + | J c i| + J c2 ] 2 - [4J afc7 (k) 

± {(| Jell + Jc2) 2 - 4J cl J c2 COS^J 1 / 2 ] } 



1/2 



•ID 



where the dispersion due to the AF exchange J a b coupling 
in ab planes depends on the two-dimensional structure 
function, 
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FIG. 6: (Color online) Spin-wave dispersions Wk (full lines) 
as obtained in the LSW theory for the parameters motivated 
by experiment^ (for J = 40 meV): (a) G-AF phase with 
Jab = 0.1425 J and J c = -0.1425 J; (b) C-AF phase with 
Jab = 0.0650J and J c = 0.0775J. Orbital excitations f2k 
(dashed lines) were obtained within the LOW theory for the 
G-AF phase (a), and for a disordered orbital state in G-AF 
phase (b). Other parameters: 77 = 0.13, V a = 0.3 J, V c = 
0.84J. High symmetry points as in Fig. [5] 



7(k) = -(cosfcj. +cosk y ). (4-12) 

The two branches of wo±(k) follow from the alternating 
FM J ci < Eq. g3D and AF J c2 > Eq. g3D exchange 
interactions along the c axis in a dimerized OVB state, 
shown in Fig. |3f a) . For the case of 77 = the spin waves 
extend up to - 1.62 J (Fig. [5J|. 

For the G-AF phase one finds the spin waves which 
depend on the (weakly anisotropic) AF exchange inter- 
actions given by Eqs. (|4.7[) and (|4.8[) . 

w G (k) = 2|(2J ah + J c ) - (2J ab 7 k + J c cosk z ) j . 

. (4 ' 13) 

For the numerical evaluation we ignored weak anisotropy 
of the magnetic exchange constants which follows from 
the spin-orbital model, and adopted the experimental 
isotropic parameters J a b = Jc = 5.7 meV, i.e., J a b — 
J c = 0.1425 J for J = 40 meV^ These parameters are 
somewhat lower than those which would result from Eqs. 
(14. 7p and (|4.8|) for the present value of J at 77 = 0.13, and 
give the width of the magnon dispersion close to 0.85 J, 
see Fig. EJa). 



Finally, we consider the C-AF phase with uniform FM 
interactions J c for which one finds the spin-wave disper- 
sion (for more details see Appendix |B|) 



1 /2 

J° ] (k) = 2{[2J Qb + I J c \(l + cosk z )] 2 - (2J afc7k ) 2 } . 

(4.14) 

This result corresponds to an idealized structure when 
the observed alternation of stronger and weaker FM in- 
teractions along the c axis (see Sec. IV Bp is ignored. 
Taking again the experimental exchange constants: 36 
J a b = 2.6 meV and J c = 3.1 meV, i.e., J a b = 0.065J 
and J c = 0.077J for J = 40 meV, one finds that the 
spin wave spectrum extends up to 0.57 J, see Fig. [fjb). 
Therefore, due to the observed strong reduction of the 
exchange interactions^ the overall width of the magnon 
band is lower in C-AF phase, while the theory predicts^ 
here a wider magnon band for rather similar exchange in- 
teractions in both C-AF and G-AF phase, as they follow 
from Eq. ([23]) . 
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C. Orbital excitations 

In contrast to spins which show particular types of long 
range order in various phases of Fig. [31 the orbitals are 
in the first instance disordered due to a robust tendency 
towards strong ID fluctuations of a and b orbitals along c 
axisj 57 ' 58 This may change, however, when lattice distor- 
tions [which induce intersite orbital interactions in 7i or b 
in Eq. (|2.18|) ] contribute and stabilize a particular type 
of AO order. Therefore, one has to employ different ap- 
proaches to determine the orbital excitations — they de- 
pend on the parameter regime and on the underlying 
orbital phase. 

First, in the OVB phase the orbitals are entirely disor- 
dered, and one has only the short range order of orbital 
singlets along c axis which imposes the dimerized mag- 
netic phase, see Fig. [3Ja). Under these circumstances, 
one finds an orbital triplet excitation for each singlet 
bond along c axis which supports local FM spin order, 

n (k) = j(ri + i« c ). (4.15) 

As these orbital excitations are local, they are disper- 
sionlcss and involve no further quantum correction to the 
energy E vb given by Eq. (|3.3|) . 

Second, although we will show below that a and b or- 
bitals are to some extent disordered in C-AF phase, weak 
long range order survives in the relevant range of param- 
eters near r) ~ 0.13, so we may start with a classical 
G-AO order at T — 0, and make an expansion around 
this state using Gaussian fluctuations. In this approach 
one rotates first the orbital operators on the sites occu- 
pied by b orbitals (down pseudospins) with (rf ) = — | by 
angle it with respect to pseudospin x axis, which leads to 
the canonical transformation: 

rf^-rf, 7?=>-7f. (4.16) 

Next we introduce a similar expansion to that considered 
above for the spin operators,"* 2 and express the orbital 
operators in terms of the respective Holstein-Primakoff 
bosons {bi,b\}, 

rr~bl 7? = i-6fo. (4.17) 

Here we assumed a robust G-type 00 (G-AO) state 
which may be used as a classical state to determine the 
orbital excitations by performing a Gaussian expansion 
around it. When only the leading terms are kept within 
the LOW theory,— one finds after the Fourier transfor- 
mation and the subsequent Bogoliubov transformation 
the orbital-wave energy, 

n { °\k z ) = J {A 2 + r\ sin 2 k z } 1/2 . (4.18) 
The spectrum has a gap at k z = 

A = { [»j(n + r 8 ) + v ] [2n + v(n + r 3 ) + «o] } 1/2 

(4.19) 



at k z =0, where 

v = 2v a - v c . (4.20) 

Note that the gap A depends on a linear combination of 
orbital interactions vq for the present form of Eq. (|2.19p , 
so the interactions along c axis and the ones in ab planes 
partly compensate each other in Eq. Q4.20[) . In fact, Eq. 
(|4.18D reproduces the earlier result obtained for vo = in 
Ref. |42j, but in general both types of orbital interactions 
originate from different distortions and are thus indepen- 
dent from each other. The orbiton dispersion demon- 
strates that the present phase is stable at finite r\ only as 
long as A > 0, i.e., in a range of vq > —r)(ri + r^). The 
orbital-wave dispersion (|4. 18|) follows from the quantum 
fluctuations along the c axis, and thus depends only on 
the zth momentum component k z . We emphasize that 
the orbital excitations are typically at higher energy than 
the spin excitations as the orbital gap A is finite, see Fig. 

EDO- 

Here we also give the values of the orbital correlations 
which enter Eqs. (|4.1|) and (I4.2[) . It is convenient to 
introduce the following integrated quantities: 




(4.21) 
(4.22) 



where A — r± + 2rj(r\ + r%)(T z ). The orbital correlations 
and the OO parameter, 

(t z ) = \-s 2 , (4.23) 

are reduced by quantum fluctuations along the c axis and 
are determined self-consistently. One finds that weak OO 
appears at finite ij > 0.08 (at r\ < 0.08 the orbitals are 
disordered and (r z ) = 0), and (t z ) ~ 0.26 for rj = 0.13, 
i.e., the OO is only about half of the classical value for the 
realistic parameters of cubic vanadates^ This demon- 
strates that the {a, b} orbitals strongly fluctuate and the 
G-AO order in rather weak. Strong orbital fluctuations 
can be also verified by calculating the intersite orbital 
correlations: 

(n-fj) = --- — [s 1 +r ] {r 1 +r 3 ) S2 ] . (4.24) 
4 r\ 

Indeed, one finds a rather low value of (ri-fj) ~ —0.428 
(not so far from the Bethe ansatz result -0.4431 for the 
AF Hcisenberg chain), and the dominating contribution 
comes not from the static term (r?r?) = — (t z ) 2 ~ 
—0.068, but from the fluctuating part, (rf tJ + rfrj) ~ 
-0.36. 

Finally, the opposite situation is found in the G-AF 
phase, where structural distortions observed below TV2 
suggest that the G-AO order sets up. In this case the 
a and b orbitals repeat each other along the chains in c 
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direction, and alternate in ab planes [Fig. [3(c)]. This 
robust C-type orbitally ordered state may be used to 
determine the orbital excitations employing the LOW 
theory^ We used again a rotation of 'down' pseudospins 
as in Eq. (|4. 1 6[) in order to obtain a uniform ferro-orbital 
state, and next expressed the orbital operators in terms of 
the Holstein-Primakoff bosons {bi, b]}, using Eqs. (|4.17[) . 
By applying a similar procedure to that used above for 
the G-AO phase, i.e., keeping only bilinear terms in the 
leading LOW order, and employing subsequent Fourier 
and Bogoliubov transformations, this leads to the orbital 
waves in the C-AO phase, with dispersion 

G (k) = J (rjn cos k z +v c + 2v a ), (4.25) 

characterized by a large gap of <~ (V c +2V a ), shown in Fig. 
UKa). We emphasize that the interactions with the lattice 
are here of crucial importance and generate a large gap, 
while the orbital gap found in the C-AO phase follows 
predominantly from the superexchange interactions and 
is therefore typically much smaller than the one in the 
G-AO phased 

D. Zero temperature phase diagram 

In order to investigate the relative stability of the mag- 
netic phases shown schematically in Fig. [3J one has to 
determine the quantum corrections due to magnetic and 
orbital excitations. The quantum corrections due to or- 
bital fluctuations were already included in the energies of 
the OVB ([33]) and G-AF ([31]) phases, where the orbital 
singlets along c axis dominate and are responsible either 
for the orbital disordered state or for weak G-AO order, 
respectively. The quantum correction to the energy of 
the G-AF phase due to the almost dispersionless orbital 
waves Eq. (|4.25[) is rather small and will be neglected 
below i£i 

The remaining quantum corrections to the classical en- 
ergy of the Neel state due to spin excitations can be 
found using the standard approach of the LSW theory. 
At T = the total energy, 

E M = E$ - 5 Em , (4.26) 
is lowered by the quantum fluctuation contribution 31 

SE M = 2J ab + | J c | - — ^ J d 3 k w M (k), (4.27) 

where label M = 0, G, G stands for a given magnetic 
phase considered here, either OVB, or G-AF, or G-AF; 
while lumO*-) in Eq. (|4.27[) is the spin wave dispersion in 
this phase. We have evaluated quantum corrections using 
Eq. P~2"7) for all three magnetic phases: OVB, G-AF 
and G-AF. It is instructive to investigate first the energy 
dependence on Hund's exchange interaction, shown in 
Fig. [7] As the quantum corrections which result from 
spin excitations are similar for all three AF phases, the 
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FIG. 7: (Color online) Energies of different phases for increas- 
ing Hund's exchange r\: C-AF phase (solid line), OVB phase 
(dashed-dotted line), and G-AF phase (long-dashed line). 
The energy of the G-AF phase is shown for 2V a = V c = 0.45 J 
(the other energies do not depend on V c ). A constant energy 
term — 2 J was neglected in all phases. The circles show for 
a comparison the energy obtained for decoupled FM chains 
along c axis, with orbital correlations described by the ID 
pseudospin Heisenberg model. 



qualitative picture obtained with these corrections and 
presented in Sec. IIII Bl is confirmed: the G-AF is stable in 
a range of realistic values of Hund's exchange r\ ~ 0.13 for 
small orbital interaction parameter V c , while increasing 
this interaction results in a transition to the G-AF phase, 
where the magnetic energy is gained on all the bonds 
after the orbitals have reoriented to the G-AO order, see 
Fig.[3J 

A transition from the OVB phase to the G-AF one un- 
der increasing rj is rather intricate. At small values of ij 
when the OVB phase is still stable, the competing phase 
with G-AF spin order is the orbital disordered phase, as 
the orbital superexchange interactions in ab planes are so 
weak (and the orbital interactions cancel each other on 
the mean field level for 2V a = V c ) that the ID pseudospin 
interaction along the c axis dominates^ the behavior of 
the orbital chain (see Fig. [7]). However, at r\ ~ 0.13 
one finds that weak G-AO order is stabilized by Ising or- 
bital interactions along the bonds in ab planes. However, 
the orbital fluctuations are still very strong in this state 
as described in Sec. IIVC1 Of course, the G-AO order 
could be further stabilized and become of more classical 
character— when 2V a > V c , but this picture of the G-AF 
phase contradicts recent experiments^ 

By comparing energies of all three magnetic phases at 
T = one finds the phase diagram of Fig. [8] To simplify 
the discussion we have adopted here the parametrization 
V c = 2V a = 2V. In fact, one expects that the GdFeOa- 
like distortions are responsible for stronger orbital inter- 
actions along c axis, and the parameter V c plays a more 
important role (than 14) in stabilizing the G-AO order 
which supports the G-AF spin order. The OVB phase is 
stable for small values of r\ and V, while for sufficiently 
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FIG. 8: Mean-field phase diagram of the spin-orbital model 
(|2.18[) as obtained for cubic vanadates in the (77, V) plane at 
T = for V c = 2V a = 2V. At the spectroscopic value of 
r\ ~ 0.13 two phases are possible: C-AF phase (for V < Vb) 
and G-AF phase (for V > V ), with V ~ 0.43 J; these two 
AF states are observed at low temperature in LaVOs (C-AF) 
and in YVO3 (G-AF), respectively. 

large V the G-AF phase takes over. At larger values of 77 
two AF phases observed in the cubic vanadates^ C-AF 
and G-AF phase, compete with each other. The range of 
stability of the G-AF phase increases with increasing 77 
as the FM interaction along c axis is then enhanced, see 
Fig. Efb). In contrast, both AF exchange interactions 
in the G-AF phase are reduced, so this phase has to be 
stabilized by larger orbital interaction V . 

V. SCENARIO FOR YVO3 
A. Peierls orbital dimerization 

Before we address the experimental situation in YVO3, 
we demonstrate an intrinsic instability of the ID spin- 
orbital chain towards dimerization. 57 In contrast to the 
ID Heisenberg antiferromagnet with fixed exchange in- 
teractions on each bond (i, the orbital interaction in 
the present case are [in the leading order, see Eq. (|2.11l) ] 
given by 

Jorb(i, i + 1) = ^(1 + 2 vn ) (Si ■ S t+1 + l) , (5.1) 

i.e., for each bond the orbital interaction is tuned by the 
spin correlation function on this bond. While at temper- 
ature T — the spins are (almost) fully polarized and 
(Si ■ Si+i) ~ 0.96*21 the spin correlations could in princi- 
ple alternate between stronger and weaker FM bonds at 
finite temperature T > 0, and then the orbital interaction 
would be modulated as follows 

J r orb(i,*+l) = ^(l±«o), (5-2) 

between even and odd bonds. Note that J Q stands here 
for the average value that will gradually decrease with 



increasing temperature. This additional temperature de- 
pendence complicates somewhat the picture of the G-AF 
phase. 

Assuming the alternating orbital interactions l|5.2p and 
performing the transformation to fermions for the cor- 
responding XY model in the orbital sector, one finds 
the following spinless fermion problem using the Jordan- 
Wigner transformation^ 

Hxy{5 ) = ± *°) (///i+i + fl+Ji)- (5-3) 

i 

The diagonalization of Hamiltonian (|5.3[) gives the energy 
spectrum of a dimerized fermionic chain, 

cos 2 k + S 2 sin fc, (5.4) 
and the total energy at T = 0: 

E(5 ) = -J — / dke„{k). (5.5) 

The energy -0.4776 J obtained from Eq. ([575)) at S D = 
is slighly lower than the Bethe ansatz result (— 0.4431J), 
while at 5 a = 1 the exact result found for the orbital 
singlets on every second bond is rigorously reproduced. 
Therefore, Eq. (|5.5|) may be considered to be a reason- 
able interpolation formula which allows one to investi- 
gate the dimerized orbital chain in the entire regime of S a . 
While an average value of the orbital correlation function 
{fi -Ti+i) increases with 5 , the chain with a constant ex- 
change interaction cannot dimerize by itself. In contrast, 
the energy E(S a ) indeed decreases when the alternation 
of the orbital interactions (|5.2[) is allowed, so the chain 
does have a tendency to dimerize (Fig. [5]). 

It is quite remarkable that already a weak anisotropy 
S a in orbital interactions is sufficient to give rather dif- 
ferent orbital correlations (ft ■ t«+i) on even/odd bonds. 
These different orbital correlations can trigger the alter- 
nation in the spin correlation functions, and in this way 
the dimerized state could be a self-consistent solution of 
the spin-orbital problem at finite temperature. We em- 
phasize that even a relatively small anisotropy S T = 0.12 
in the orbital correlations, 

S T = \{n- - (r i+1 • T i+2 )\, (5.6) 

is already sufficient to generate considerable anisotropy 
in the magnetic exchange constants J c \ and J c i along c 
axis (Fig. fTUj) . The exchange constants of Fig. [TTJl were 
obtained with J = 30 meV — this reduction of the en- 
ergy scale by a semiempirical factor of 0.75 from that 
given by the analysis of the optical spectrum^ was nec- 
essary as otherwise the model (|2.18p would predict too 
large exchange constants for the G-AF phase. Further- 
more, we note that the above anisotropy S T is obtained 
already with S = 0.017 when the mapping to the fermion 
problem (|5.3[) is used (Fig. [5]). Of course, this problem 
requires a self-consistent solution at finite temperature 
as we discuss in Sec. IV Bl 
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FIG. 9: (Color online) Orbital correlation functions (fj ■ fi+i) 
at even and odd bonds along a dimerized ID chain, as ob- 
tained within the XY model for increasing anisotropy S a in 
the exchange constants, see Eq. (|5.2[l . The energy E(5 ) 
(dashed line), obtained using spinless fermions Eq. (|5.5jl . de- 
creases with increasing S , while the average energy in an 
orbital chain with the same exchange interaction J at each 
bond would increase (dotted line). 
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FIG. 10: (Color online) Exchange constants J a b (dashed line), 
Jci and J C 2 (solid lines), all in meV, as obtained for the ide- 
alized C-AF phase with condensed xy orbitals (n c = 1), Eqs. 
(|2.8[) . and orbital disordered state along c axis with dimer- 
ized orbital correlations. Vertical line indicates the value 
of T) = 0.13 estimated from the atomic data££ and from 
the optical dataiS for LaVOa. Parameters: J = 30 meV, 
(n ■ f i+ i) = -0.4431, S T = 0.12 (EH). 



B. Reduction of exchange constants by orbital 
fluctuations 



Although the value of J ~ 40 meV deduced from the 
neutron scattering data^ for YVO3 gives a consistent 
description of the temperature dependence of the optical 
spectral weight for the high-spin excitations along c axis 
in LaV03, there is a fundamental problem concerning the 
size of magnetic exchange constants, particularly in the 
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FIG. 11: (Color online) Reduction of exchange constants J a ^ 
(dashed line), J c \ and J7c2 (solid lines), all in meV, in the 
dimerized C-AF phase due to orbital bond fluctuations be- 
tween FM and AF bonds, as given by Eqs. (15.7 [) - (15.9|) . Or- 
bital disordered state along c axis is assumed at p — 0. Ex- 
perimental values of exchange constants found for the C-AF 
phase of YVO3,— shown by circle (J a t — 2.6 meV) and by 
diamonds (J c i ~ 4.2 meV and J C 2 — 2.0 meV), are nearly 
reproduced for moderate fluctuations with p = 0.30 (vertical 
dashed line). Parameters: J = 30 meV, r\ — 0.13, S s = 0.35. 



exotic C-AF phase of YVO3 , stable in the intermediate 
temperature range TV 2 <T < Tjyi. First of all, the cal- 
culations performed using the mean-field approach and 
assuming rigid OO (see Fig. [3]), as in Ref. [9|, predict too 
large exchange constants in both phases when J = 40 
meV is assumed. In fact, for the G-AF phase one finds 
then the values of both J c and J a b being larger by ~ 25% 
than the respective experimental values of Ulrich et al.£^ 
Moreover, in experiment one finds an (almost) isotropic 
G-AF phase with J c = J a b, while the present model pre- 
dicts (except at small rj < 0.10) an anisotropy between 
c axis and ab planes, with J c > J a b, see Eqs. (|4.7p and 
(|4.8[) . This suggests that already for the G-AF phase 
some 'dynamical' reduction mechanism of the magnetic 
exchange constants is at work, which we simulate by re- 
ducing the superexchange energy scale down to J ~ 30 
meV. Indeed, taking an average value of the magnetic 
exchange constants over three cubic directions we arrive 
then at the experimental result J c — J a b ^5.7 meV. 

While the above procedure could be still considered 
as a fair agreement between the theoretical model and 
experiment, it is surprising that the magnetic exchange 
constants in the G-AF phase cannot be obtained from 
the model using the same parameters. In fact, the values 
of J c and J a b shown in Fig. [TJJ] for r\ = 0.13 are by 
almost a factor of 2 larger than those deduced from the 
neutron scattering data at 85 K.— This strongly suggests 
that some of the assumptions used so far to derive the 
values of and J% b from Eqs. I|4.ip and P~2"|) have to 
be reconsidered. 

One of the most puzzling experimental features in 
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YVO3 is the nature of the structural transition at T s , 
which removes the orbital degeneracy and induces the 
splitting A between the xy orbitals and the yz/zx dou- 
blet, see Fig. [1] We anticipate that this splitting is not 
large enough to impose strict freezing of charge in xy 
orbitals. Thus we expect that some orbital fluctuations 
should still be present in the intermediate temperature 
regime T^ 2 < T < T s before the orbitals undergo the 
transition into the C-AO phase (supporting G-AF spin 
order) below Tjv2, as shown in Fig. E^c). 

Qualitatively, we illustrate the consequences of orbital 
fluctuations on the magnetic exchange constants by con- 
sidering a plaquette which includes two bonds along c 
axis and two other bonds along either a or b axis. If the 
c orbitals are occupied at each site, and a/b orbitals fluc- 
tuate, a representative state of such a plaquette contains 
4 electrons in c orbitals, and 2 in each of two other states, 
a and b. The effective superexchange Hamiltonian (12. 9p 
contains the terms with double orbital excitations on the 
bonds, cx T~iTj, see Eq. (|2.16|) . Such terms on the bonds 
along a (or b) axis generate a/b orbital configurations on 
each site i and j. Only one of these two orbitals (a or 
b) is active along this particular bond, and it resembles 
the bond along c axis before the orbital fluctuation took 
place. As a result, such fluctuations lead to (locally) FM 
contributions in the ab planes, and to (locally) AF con- 
tributions along c axis — both of them will reduce the 
actual values of and J^ b exchange constants. 

Following the above idea, we introduce effective mag- 
netic exchange constants, 
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Jab 



(1 + 5.) [(1 - p) J c c (0) + pjg(0)] , (5.7) 
(1 - S s ) [(1 - p) J c °(0) + pjg(0)] , (5.8) 
(l-p)J Q c b (0)+pJ c c (0), (5.9) 



as a superposition of two contributions obtained for the 
undimerized state without xy orbital fluctuations (for 
n c = 1), J,F(0) and J^(0), calculated as described in 
Sec. IIV Al The probabilities (1 — p) and p refer to the 
initial state with c orbitals occupied (n c = 1), and to 
the configuration with flipped orbitals after the plaque- 
tte fluctuation has occurred (n c = 0), respectively. The 
result of the numerical calculation for the usual parame- 
ters shows that one arrives almost at experimental values 
of the magnetic exchange constants when moderate or- 
bital fluctuations with p — 0.30 considerably reduce the 
exchange constants (see Fig. [TTj) . For the experimental 
anisotropy 6 S one finds large alternation of the FM ex- 
change constants along c axis with respect to the average 
value, 



# = (l-p)#(0)+p# 6 (0). 



(5.10) 



Next, we analyze the spin excitations in the dimcr- 
ized C-AF phase in order to calculate the spin correla- 
tions, the quantum fluctuation correction to the ground 
state energy (see Appendix [B} , as well as the free energy 
at finite temperature, see Sec. IV CI The effective spin 



25 



20 



> 


15 j- 


CD 








10 ' : 


3 






5 - 




FIG. 12: (Color online) Spin-wave dispersions Uk (full lines) 
as obtained in the LSW theory along the representative di- 
rections in the Brillouin zone for the dimerized C-AF phase 
with experimental exchange constants; 3 ^ J a ^ = 2.6 meV, 
J c = 3.1(1 ± S a ) meV and S a — 0.35. These interactions 
are obtained by considering plaquette fluctuations of spin ex- 
change interactions as described in the text (see also Fig. 1111) . 
Parameters: J = 30 meV, i) = 0.13, 8 S = 0.35, p = 0.30. The 
experimental points of Ref. [3^ measured by neutron scat- 
tering at T — 85 K are reproduced by circles (the effective 
linewidths are not shown). The high symmetry points are: 
T = (0,0,0), M = (7r,7r,0), R= (7T,7r,7r), Z = (0,0, tt). 



Hamiltonian for this phase is given as follows: 
Us = J?(1 + 6 S ) S 2l -S 2l+1 

(2i,2i+l)||c 

+ Jf(l-<5 S ) S 2i -i-S 2 i 

(2i-l,2i)||c 

(ij)\\ab 



(5.11) 



Following the LSW theory, the spin wave dispersion is 
given by 



wo±(k) = 2{ [2J ab + \J C \ 

± ^(cos^ + ^sin 2 ^) 1 / 2 ] 2 - (2x67k) 2 } 



1/2 



(5.12) 

For the numerical evaluation of Fig. [12] we have used the 
experimental exchange interactions^ J a b = 2.6 meV, 
J c = 3.1 meV, S s = 0.35. Indeed, large gap is found 
between two modes halfway in between the M and R 
points, and between the Z and T points (not shown), 
respectively. Two modes measured^ and obtained from 
the present theory in the unfolded Brillouin zone follow 
from the dimerized magnetic structure. 

The microscopic reason of the anisotropy in the ex- 
change constants J c \ and J c2 is the tendency of the or- 
bital chain to dimerize, as we have demonstrated in Sec. 
IV Al Such a dimerized orbital chain may only be stable, 
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FIG. 13: (Color online) Result of the self-consistent calcu- 
lation of intersite correlations in the dimerized C-AF phase 
along c axis for increasing temperature: (a) orbital (fj ■ ri+i), 
and (b) spin (Si-Si+i). The correlations on stronger (weaker) 
FM bonds are shown by solid (dashed) lines. Dashed-dotted 
line in (b) shows the order parameter (S z ) in the C-AF phase 
as obtained from the LSW theory. Parameters: n — 0.13, 
J = 30 meV, and p = 0.30, see Fig. ITD 



however, if the corresponding interactions in the orbital 
sector (|5.2|) alternate, i.e., 5 > 0. This becomes possible 
at finite temperature when also intersite spin correlations 
may alternate along the c axis, supporting such a dimer- 
ized state. Although a completely satisfactory treatment 
of the spin correlations in a broad temperature regime 
which covers the symmetry broaken C-AF phase is not 
possible at the moment, we have employed the LSW the- 
ory to calculate the spin correlations (Si ■ Sj+i), as ex- 
plained in the Appendix [Bj 

The result of the self-consistent calculation of spin and 
orbital correlations along the c axis in the dimerized C- 
AF phase is shown in Fig. [T3J The driving force to 
stabilize the dimerized state is the instability of the or- 
bital chain which leads to rather strong anisotropy in the 
orbital correlations [Fig. [T3Ta)]. On the contrary, the 
spin correlations differ only by a rather small amount 
(unlike in the OVB phase), as the large spins 5 = 1 are 
far less susceptible to follow the dimerized structure, and 
the long range spin order is supported by the exchange 
interactions in all three directions. The energy of the 
dimerized state is lower than that of the undimerized C- 
AF structure. Apparently, a weak anisotropy between 



(Si ■ S i+ i) ~ 0.93 and ~ 0.79 on stronger/ weaker FM 
bonds encountered at T — 77 K [see Fig. [Tffib)] is already 
sufficient to trigger a phase transition to this phase from 
the G-AF phase stable below T^2- Why this transition 
may really happen in YVO3 we explain in the following 
Section. 



C. Mechanism of the phase transition from G-AF 
to C-AF phase 

The transition from G-AF to G-AF phase in YVO3 is 
puzzling as the magnetic order changes completely at fi- 
nite temperature Tjv2 — 77 K, and the magnetic moments 
reorient^ The observed change of the spin and orbital 
pattern indicates that the spin-orbital superexchange in- 
teractions are frustrated and it is easy to tip the balance 
of these interactions and to change completely both the 
magnetic and orbital order. As the transition between 
the two phases occurs at finite temperature, the entropy 
has to play an important role, so it was suggested before 
that the large orbital entropy due to orbital fluctuations 
in the G-AF phase could be released at Tjv2 and trigger 
the transition^ A closer inspection of the present model 
and the reconsideration of recent experiments show, how- 
ever, that the situation is somewhat more intricate. 

First of all, we have already emphasized that the mag- 
netic exchange constants are reduced in the G-AF phase, 
and we presented a possible mechanism responsible for 
this reduction in Sec. IV Bl As a result of orbital fluc- 
tuations, the average energy of magnetic excitations is 
lowered in the G-AF phase (Fig. ^ , so one expects that 
the spin entropy might play an important role as well. 
Using the spin and orbital excitations derived already for 
both phases in the previous Sections, we estimate these 
entropy contributions assuming that the excitations are 
independent from each other. The spin waves are given 
by Eqs. (|4.13p and (|5.12[) . while the orbital excitations 
by Eqs. (|4.25[) and (|4.18p . Here we will ignore the change 
of the orbital excitations in the dimerized G-AF phase as 
this gives only a marginal contribution to the entropy 
of the G-AF phase, and does not influence the magnetic 
transition at TV2 significantly. 

The spin and orbital entropy normalized per one vana- 
dium ion is calculated using standard formulae: 



S c 



**^£{]og(l- e HW0; 

k 

log (1 _ e -0«c-eoU 

fe^i-5>g(l-e-^«), (5.13) 



s o = fc B ri£{i g(i-e-^ k )) 



+ k B T 



1 _ e -pn a (k) 



(5.14) 
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FIG. 14: (Color online) Entropy of the C-AF (Sc) and G -AF 
(Sg) phase as obtained for the spin-orbital model (|2,18[l us- 
ing the experimental values of magnetic exchange constants 
in both phases. The dominating contributions result from 
spin excitations (dashed-dotted lines), while the orbital con- 
tributions (dashed lines) are much smaller, but give also a 
higher entropy in the C-AF phase. Parameters: J = 40 meV, 
■q = 0.13, V a = 0.30J, V c = 0.84J. 



where (3 = 1/ksT, and N (Ni) is the number of k (k) 
values. The entropy consists of the spin and orbital en- 
tropy terms for each phase. All summations are over the 
Brillouin zone which corresponds to the undimcrizcd C- 
AF phase. Using the parameters consistent with the ex- 
perimental data of Ulrich et al££- one finds (see Fig. [M]) 
that: (i) the entropy Sc for the C-AF phase is larger that 
Sg for the G-AF phase, and (ii) the spin entropy grows 
significantly faster with temperature than the orbital en- 
tropy for each phase. Therefore, we conclude that the 
spin entropy gives here a more important contribution 
and decreases the difference between the free energies of 
both magnetic phases in the temperature range T ~ Jjv2- 
It has been argued before^^ that the difference be- 
tween the energies of both phases, Eq and Ec, has to be 
small at T = 0. Indeed, we evaluated the free energy of 
both phases using the above entropies (|5.13|) and 1)5 .14p . 

T c = Ec-TSc, (5.15) 
•F g = E G ~TS G , (5.16) 

and found that Ec — Ec — 1 meV, and the transition 
from G-AF to G-AF phase is reproduced at the exper- 
imental value of the temperature Tjv2 when the orbital 
interactions are chosen properly. In Fig. [15] we show a 
representative case with V a = 0.30J, V c = 0.84J with 
J = 40 meV. Of course, this fit is not unique and V a 
(V c ) could be somewhat smaller (larger), but the energy 
difference Ec — Ec at T = remains close to 1 meV in 
all cases. Note, however, that too large values of V c are 
not allowed, as then the G-AF phase gets destabilized by 
orbital excitations (|4.18p . 



-19 




50 100 150 



T(K) 

FIG. 15: (Color online) Free energies of the G-AF (Tc, solid 
line) and G-AF (Tg, dashed line) phase as obtained for the 
spin-orbital model (|2.18[1 using the experimental values of 
magnetic exchange constants^ in both phases. The exper- 
imental magnetic transition temperatures, Tjv2 — 77 K and 
Tjvi — 116 K, are indicated by arrows. Parameters are the 
same as in Fig. 1141 



VI. SUMMARY AND CONCLUSIONS 

The present study shows the importance of ti g orbital 
degrees of freedom in cubic vanadates. We have pre- 
sented the spin-orbital model for cubic vanadates and an- 
alyzed its possible solutions in various parameter regimes, 
using extensively the decoupling of spin and orbital de- 
grees of freedom. Although the model is more general, 
we have focused on the solutions which arise in the case 
of anisotropic occupancy of ti g orbitals, with xy orbitals 
singly occupied at each site. This state is believed to be 
realized in cubic vanadates, at least in low temperature 
phases with magnetic long range order. When Hund's ex- 
change and orbital interactions promoted by the lattice 
are weak, the superexchange is strongly frustrated and 
gives a rather exotic dimerized OVB state, with orbital 
singlets alternating along the c axis and stabilized at ev- 
ery second bond by ferromagnetic spin correlations. In 
this way, spin and orbital correlations support each other 
and demonstrate a unique instability of the spin-orbital 
system towards a dimerized statei 37 ' 57 This instability 
turns out to play also an important role at finite temper- 
ature in YVO3, but in a different regime of parameters 
where its mechanism is more subtle. 

When Hund's exchange or the orbital interactions in- 
crease, the OVB ground state is disfavored and a par- 
ticular type of long range magnetic order emerges in- 
stead from the frustrated superexchange interactions in 
cubic vanadates. These other AF states (G-AF and G- 
AF phases), as well as the OVB state itself (at low Jh), 
demonstrate a close interrelation between magnetic and 
orbital order, with complementary behavior of spin and 
orbital correlations, known as the so-called Goodenough- 
Kanamori rules^jM While in some cases these rules (and 
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the underlying decoupling of spin and orbital operators) 
work well, we have presented the case of the G-AF phase 
with rather disordered orbitals, where it is likely that 
joint spin-orbital fluctuations also play a role^ and it 
would be necessary to include them for a more quantita- 
tive comparison with experiment. 

A detailed analysis of the possible solutions of the spin- 
orbital superexchange model supplemented by the orbital 
interactions induced by the lattice demonstrates that two 
different types of AF order, C-AF and G-AF phase, com- 
pete with each other in the parameter regime relevant for 
YVO3. However, the energetic proximity of these two 
phases in a particular parameter regime could explain 
possible changes of magnetic order by pressure or mag- 
netic field — when the microscopic parameters are fixed, 
one or the other phase could be stable at T — > 0. The 
situation changes at finite temperature, however, when 
the spin and orbital excitations are of importance and 
may tip the energy balance between given two types of 
order by the entropy term. In fact, we have shown that 
this is likely to be the microscopic explanation of the ob- 
served first-order phase transition and switching of the 
magnetic order in YVO3 at Tjv2- 

Our study has established that the nature of the tran- 
sition from the G-AF to G-AF phase at T/v 2 observed in 
YVO3 is complex and several factors have to come to- 
gether to trigger it when temperature increases: (i) the 
presence of active t2 g orbital degrees of freedom opens a 
possibility of two different types of AF order which may 
compete with each other; (ii) rigidity of the G-AO order 
in the G-AF phase hampers possible free energy gains 
when spin or orbital excitations are created (as spin in- 
teractions are rather strong and the orbital gap is quite 
large); (iii) the change of structure observed at Tn2 not 
only helps to stabilize the weak G-AO order, but also re- 
leases more orbital fluctuations when the xy orbitals be- 
come active and their occupancy is not fixed — such fluc- 
tuations result in turn in fluctuating magnetic exchange 
constants and lead to the reduction of the characteris- 
tic energy scale for the spin excitations, and finally (iv) 
the spin correlations have to be weakened by increasing 
temperature to participate in a joint spin-orbital dimer- 
ization in the G-AF phase. Thus, the difference between 
the G-AF and G-AF phase of YVO3 is much deeper than 
simply the observed difference in the magnetic order. It 
is far more important that the orbital state softens at the 
transition at Tjv2 to the G-AF phase and this change hap- 
pens in a concerted way with the observed reorientation 
of the magnetic moments. In addition, the intrinsic insta- 
bility in the orbital sector towards dimerization, which is 
incompatible with the magnetic order in the G-AF phase 
and is blocked by spin correlations in the G-AF phase 
at T = 0, becomes possible when the intersite ferromag- 
netic spin correlations along c axis have been somewhat 
weakened with increasing temperature. 

Although we have suggested a plausible scenario of the 
observed magnetic phase transition in YVO3, the micro- 
scopic theory of the G-AF phase at finite temperature 



remains still to be constructed. It is not clear at the 
moment to what extent the xy orbital fluctuations are 
released in the intermediate magnetic phase and are still 
present up to the structural transition at T s ~ 200 K (see 
Fig. 1). It could well be that spin-orbital entanglement 
in excited states plays a role in this temperature regime 
and prevents reliable evaluation of the magnetic exchange 
constants in the G-AF phase using the conventional de- 
coupling of spin and orbital operators. Furthermore, it is 
puzzling whether dimerization also plays a role in reduc- 
ing the the magnetic order parameter in the G-AF phase 
which is hard to explain using the spin-wave theory, or 
the above entanglement is the main reason responsible 
for this reduction. Note however that it could be argued 
that the observed orientation of the magnetic moments 
which are close to lying within the ab planes is enforced 
by the dimerization in the G-AF phase. 

Some other problems remain still open and should be 
treated in future more complete theory. We note that 
also in the G-AF phase a considerable reduction of the 
magnetic order parameter— goes beyond that expected 
from the quantum fluctuations^ We believe that the 
relativistic spin-orbit coupling cx XLiSi contributes sig- 
nificantly to the magnetic properties in the entire regime 
of temperature^ in particular also to the spin correla- 
tions in the G-AF phase, and could reduce the observed 
value of the magnetization. In fact, it would also break 
the symmetry in the spin space and determine an easy 
axis for the AF order parameter. In the present study the 
spin-orbit coupling A was ignored, as in the considered 
regime of J ^> A it could lead just to the perturbative cor- 
rections of the presented spin and orbital excitations. In 
contrast, in the regime of A ~ J it would lead to ordering 
of orbitals with complex coefficients, {\xz) ± i\yz))/>/2, 
with finite orbital angular momentum^ Although the 
cubic vanadates are not in this regime of parameters, 
finite spin-orbit coupling A would be crucial for quanti- 
tative understanding of spin (and orbital) excitations in 
the entire parameter regime. This interaction provides 
another mechanism for the softening of spin excitations 
in the G-AF phase, which would complement the scenario 
considered in this paper. We also note that a small G- 
like magnetization component was observed as well in the 
G-AF phase in the temperature range T/y 2 < T < Tjvi. 
Therefore, it is likely that this magnetic phase is still 
more complex than suggested in the present paper, and 
requires a more careful analysis. Recent progress in ex- 
perimental methods makes it possible to measure also 
orbital excitations ) 63 ' 64 and information on the orbital 
excitations in YVO3 would be instrumental to resolve 
some of the above problems. 

Summarizing, we have presented the consequences 
of the microscopic spin-orbital model in the parameter 
regime relevant for cubic vanadates and suggested a sce- 
nario which explains the magnetic transition between the 
G-AF and dimerized G-AF phases observed in YVO3. 
This study indicates a close relationship between the ob- 
served magnetic correlations in the ground state and the 
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structural transition, which in case of YVO3 occurs well 
above the first magnetic transition. Thus we conclude 
that a careful analysis of the mechanism of the struc- 
tural transition and its dependence on the actual chem- 
ical composition is challenging and needed for complete 
theoretical understanding of the experimental phase dia- 
gram of cubic vanadates. 
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APPENDIX A: DERIVATION OF THE 
SPIN-ORBITAL MODEL 

The effective superexchange interactions between two 
V 3+ ions in d 2 configuration with spin 5 = 1 (triplet 3 T 2 
state) at sites i and j for a bond (ij) oriented along one 
of the cubic axes 7 = a, 6, c originate from virtual charge 
excitations by the hopping processes which involve two 
active t 2g orbitals along its direction. As an example, we 
consider here a bond along c axis (7 = 0), with active a 
(yz) and b (xz) orbitals. In this case the charge excita- 
tion by either a or & electron leads to one of three possible 
d 3 excited states: abc, a 2 c or b 2 c (see Fig. [5]). The actual 
configuration c 1 of the inactive orbital c enters via the 
constraint (|2.6[) . The total spin per two sites is conserved 
in the dfd 2 — > dfdj excitation process, i.e., the electron 
transferred in the excitation process and two other elec- 
trons on the d 3 site are either in high-spin (5=1) state, 
or in low-spin (5 = 5) state. Therefore, when the sec- 



ond order processes dfd 2 



dH) 



dfd 2 are analyzed, 



one has to project the d 3 configuration generated after 
an individual hopping process on the respective d 3 eigen- 
states. Similar, when a deexcitation process took place, 
one has to project the resulting df configuration on the 
triplet 3 T2 ground state. 

The general form of the effective Hamiltonian follows 
from symmetry considerations for the possible dfd 2 — > 
dfd]j — > dfd 2 processes which contribute to the superex- 
change. The total spin states in the excited states are 



well described by the spin operators: 

Pus (4 S,) = s , 2 g + 1) & ■ % + 2), (Al) 
Pls(4 Sj) = s{2 s + 1) (Si • S, - 1), (A2) 

which correspond to the high-spin Vns(Si, Sj) and low- 
spin Vhs{Si, Sj) excited states, respectively. The orbital 
state is described by the orbital operators Q n (i,j), where 
n refers to different excited states*^ In the present case 
n = 1 corresponds to high-spin 5 = | excited states in 
Fig. [2{a), 71 = 2 to low-spin S — \ excited states in Fig. 
[2jb), while n = 3 describes the orbital state realized for 
the excitations of double occupancies shown in Fig. [DJc). 
As the orbital quantum number is conserved along the 
hopping process (|2.2[) , either the same two orbitals are 
occupied before and after the virtual excitation, or an 
orbital fluctuation shown in Fig. [2] takes place, and the 
occupied orbitals are interchanged between sites i and j. 
The latter processes are unique for the t 2g orbitals and do 
not occur for degenerate and singly occupied e g orbitals, 
where the orbital quantum number is not conserved and 
single orbital excitations are possible instead^ 

In the case of cubic vanadates one arrives therefore at 
a general expression, 



"=£{- 

(ij) 

1 t 2 



1 i 2 



: (S i -S j + 2)Q 1 (i,j) 



3e( 4 A 2 )' 

(Si • Sj-l)Q 2 (i,j) 

(£ • 5, -1)23(^)1, (A3) 



3e( 2 E) 

1 / t 2 t 2 



2 VeTO e( 2 T 2 ) 



The first term cx t 2 jei^A-i) is FM, while the remaining 
terms stand for different AF contributions. The coeffi- 
cient 1/2 in the contributions due to 2 T\ and 2 T 2 excited 
states follows from the projection of the double occupan- 
cies of one of the active orbitals, either afci or bfci [Fig. 
|2fc)], onto the eigenstates of V 2+ ions. The orbital states 
which contribute to the above structure of superexchange 
(|A3j) depend on the bond direction; here we give as an 
example a complete expression for the bonds (ij) along 
the c-direction. 
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1 t 2 



1 2t 2 



3 e(±A 2 ) 



(Si ■ Sj + 2) [(1 - n io )(l - n jfc ) + (1 - n lb )(l - n,- a ) - (albMaj + 6ta j at6 i )n ic n ic ] 



3 e( 4 A 2 ) 
1 2i 2 



(5.; • Sj + 2) [(1 - n ic )n jc + n ic (l - n JC )] 



1 



i e 



3 e( 2 E) 



(Si ■ S 3 - 1) [(1 - n lc )n jc + n ic (l 



+ 



2 Ve( 2 Ti) e( 2 T 2 ) 



(Si ■ Sj - 1) [(1 - n ia )(l - n jb ) + (1 - »j fe )(l - Uja) ~ n jc - n jc + 2] 



2 Ve( 2 Ti) e( 2 T 2 ) 



(Si-Sj-i) (.rhn h, ■ v<X«,) }. 



(A4) 



The operators of and b\ are (spinless) fermion creation 
operators in the active orbitals \a) and \b) at site i, respec- 
tively, while n i7 = jj^fi are fermion number operators in 
state 1 7) at site i, with 7 = a, 6, c. 

The effective interactions on the bonds within the ab 
planes may be now obtained by rotating Eq. (|A4[) to the 
bonds oriented along either a or & axis. Note that the or- 
bital operators which correspond to the active \a) and \b) 
orbitals are then replaced by either \b) and |c) (for a axis), 
or by I a) and |c) (for b axis), while the general structure of 
the superexchange Hamiltonian (|A4[) remains the same. 

As both FM and AF terms are present in Eq. (|A4I) . 
the superexchange interactions are frustrated. It is in- 
structive to consider the limit of Jh — > in which the 
multiplet structure of V 2+ ions collapses to a single ex- 
citation energy U . In this case the interactions simplify 
considerably, and the terms cx n JA , which originate from 
the excitations with three different orbitals occupied at 
the same site cancel each other. This feature is analo- 
gous to the similar compensation of the high- and low- 
spin processes in the superexchange models with degen- 
erate e g orbitalsJ However, as a new feature one finds 
a nonvanishing contribution due to the orbital fluctua- 
tions, oc (aj&^Oj + H.c), as the terms which originate 
from the high- and low-spin processes add to each other. 
As usual, the double occupancies in the excited states 
lead to the AF terms as a consequence of the Pauli prin- 
ciple. These simplifications lead to the following form of 
the effective Hamiltonian along the c axis in the limit of 
Jh ->0, 

H c (t) = 0) = J % + X ) 

x (n ia rij b + riibUja + </,' //,'/« , + , 

(A5) 



expression may be also written in a more compact form, 
H c (r, = 0) = ~ J J] (Si ■ S, + Tj + \rnn ) , (A6) 

fa') He 

where the scalar product of orbital pseudospin operators 
n = > T r> T /} is defined by: 



(A7) 



Here we use spinless fermion operators a\ and b\ , but one 
could also introduce instead Schwinger boson operators. 
For the bonds along either a or b axis similar expressions 
obtained from Eqs. (|A7[) by cyclic permutations of the 
orbitals {a, 6, c} have to be used. If in addition the c 
orbitals are condensed (rii c = 1), as in YVO3, one finds 
a simplified form of Eq. (|A6|) for the bonds along c axis, 



(A8) 



<*i>ll c 



The orbital interactions are then purely classical on the 



bonds in ab plane as (tJ • fj + jriinj)"" 
bonds. 



6—1 



for these 



APPENDIX B: SPIN AND ORBITAL 
EXCITATIONS IN THE DIMERIZED C-AF 
PHASE 

Here we explain the full algebraic structure of the spin 
and orbital wave problem in the dimerized C-AF phase. 
Its solution gives the both types of excitation energies, 
and provides a systematic method to evaluate both the 
value of the order parameter, 



where J = 4t 2 /U is the superexchange interaction. This 



(s; 



s-ss z , 



(Bl) 
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and the intersite spin correlations along c axis, 

0i ■ S i+1 ) = Cu+i = Co + 5Ce 1 ^ , (B2) 

where Zi is the zth coordinate of the vactor Ri corre- 
sponding to site i. If exchange interactions alternate 
along c axis, as given by Eqs. (|5.7|) and (15. 8|) . the alter- 
nating part 8C of the intersite spin correlation function 
is finite. 

In order to evaluate the order parametr (|B1[) and the 
intersite spin correlations (|B2[) in the dimerized C-AF 
phase within the LSW formalism it is convenient to ar- 
rive first at the boson representation of the spin Hamil- 
tonian. Therefore, we performed the transformation to a 
ferromagnet ()4.9j) and the subsequent Holstein-Primakoff 
transformation (|4.10[) to {ai,a\} bosons. One finds the 
following form of the above averages, 

5S* = S-{a\a t ), (B3) 
C M+ i = S 2 -2SS Z + (atat +1 ), (B4) 

and for the quadratic (LSW) Hamiltonian 

#LSW = Jab ^2 ( Hi + n i + + a » j) 

(ij)\\ab 

+ ^ E + 

(i,t+l)||c 

x (m + n l+1 - a\ +l a t - a\a t+1 ), (B5) 

where J~ is the average value (|5.10|) of the FM exchange 
interaction along c axis. 

Next we employ the Fourier transformation to boson 
operators in reciprocal (momentum) space 

4 = -LVe ikRi aj, a lL = ^=T,e- <kat a i , (B6) 



which gives the LSW Hamiltonian in reciprocal space, 
#lsw - J2{ AJ ab a W + 2JS7(k)(4°tk + «- k a k ) 

k 

+ 2jC[(l - cosfc 2 )aj c a k + i6 s sin k z a^a k+ o\ j , 

(B7) 

where Q = (0, 0, ir) is the wave vector which corresponds 
to the doubling of the unit cell along c axis due to the 
dimerized C-AF spin structure. 

In order to find both the energies of spin wave excita- 
tions and the average values of the correlation functions 
at finite temperature T, we introduce here temperature 
Green's functions for boson operators in the momentum 
space using the notation of Zubarev^^ The first of 
them satisfies the following equation of motion, 

w«Okk»« = ^ + ((K,^lsw]|4»- (B8) 
It depends on energy w and generates three more 
Green functions: ((a k +Ql a k))^> (( a -kl a k;))"> and 
(( a -k+Ql a ic))"- I n 1S nex t convenient to introduce the 
following expressions which define the algebraic structure 
of the spin wave problem: 

A k ± = 2J c c (l±cosfc z )+4J' ( g, (B9) 
Bk = 4J^ 7 (k), (BIO) 
A k = 2J°8 s smk z . (Bll) 

The respective system of equations of motion generated 
by Eq. jB8} is: 



- w (k) iA k B k 

-iA k A k+ - u c (k) 

-B k -A k _-c^ c (k) 

-B k iA k 





-»A k 
-A k+ - wc(k) 



/ (Kl4»- \ 

((a k +q|a k »^ 

v ««-k+qi4»« y 



-1° 

2tt ( 





(B12) 



Eq. (|B12p has a typical structure obtained for elementary excitations in the random phase approximation (or in 
the LSW theory) for an antiferromagnet. One finds two positive eigenvalues o>c±(k) given by Eq. (I5.12p . and two 
negative ones, — uc±(k). 
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By solving the system of Eqs. (|B8|) one finds the following Green's functions: 

I t \\ . 1 (^ 2 - Aj + + Bp(u + A k _) - A 2 ^ - A k+ ) 
«a k |a k » w = +- {w2 _^ +(k)}{w2 _^_ (k)} > 

i t u _ * A k [(^ + A k+ )(^ + A k _) + ^-A k ] 
«« k+Q |a k )). - -- K _^ +(k)}{w2 _^_ (k)} - 



// t I t\\ 1 j? k (^ 2 -^ + + ^-A^) 



// t i t\\ , * A k £ k (,4 k+ + A k _) 

\( a -k+o a k/)^ = + 7^"T^ 5 — TTTTF 7 2 — 77TT- ( B16 ) 

k+Qi k// 2^{w 2 -^ + (k)}{w 2 -^_(k)} 



They contain complete information about the bosonic 
correlation functions which appear in Eqs. ()B3|) and 
(|B4|) . They are obtained from the temperature Green's 
functions using the fluctuation-dissipation theorem^ 



SS Z = 



sc 



k 

f L 2S(( at |4)) M -j 

^ k ^ 

^E( fl k«: 

k 



(B17) 



k+Q/ 



2»«« t _ k l4»- , 



1 



(B18) 



where (3 = l/k B T. 

In a similar way one may find the orbital excitations 
and the respective Green's functions needed to determine 
the alternation of the orbital correlations in the dimerized 
structure, 



(n ■ t 1+ i) = 7m+i = T^li + STe*"*, (B19) 
and the renormalized value of the order parameter, 



(B20) 



As in case of spin operators, we used the rotation (|4.16|) 
of orbital operators to the ferro orbital state, followed 
by the Holstein-Primakoff transformation (|4. 17|) to the 
respective boson operators {b^b\}. One finds the LOW 
Hamiltonian, 



where pi = b\b i . 

In spite of the ID nature of orbital dispersion (|4.18l) . 
the Fourier transformation to boson operators in the 
reciprocal (momentum) space is three-dimensional and 
takes here the form 



% h 



i 



(B22) 



Simalr to spin case, we find the energies of orbital wave 
excitations and the average values of the boson opera- 
tors at finite temperature T, using temperature Green's 
functions for boson operators in momentum space i 65 ! 66 
The system of equations is generated by the following 
equation of motion, 



"(K\bi)) u = ±- + ((K,H h0 w}\bt)) 
Zir 



k //w > 



(B23) 



and equations for three other Green functions: 
(K +Q K))^ {(b-M) u , and ((b 

-k+Ql 6 k))"> follow. 

The following expressions define the algebraic structure 
of the orbital problem: 



A k = 
+ 

5 k = 

e k - 



r\ + r\(r\ + r 3 ) - ~(1 - C )[n - r)(ri + r 3 ) 



-(1 + 2 vn - W3 )Y ab , 
ri-i(l-Cb)ri(l-r 
-(1 - C )S o n(l - rj) smk z 



(B24) 
(B25) 
(B26) 



#low = Jr)(n +r 3 ) ^ (p 

(ij)\\ab 
1 



Pj) 



(i,i+l)\\c 



(Pi +Pi+i) 



b\b i+1 ) 



IjCo^nCl-r?) £ i h \ + A + b\K +1 ). (B21) 

(i,i+l)\\c 



Here we used the short hand notation for the spin corre- 
lation function in ab planes, 



Y ab = (Si ■ Sj) + 1, 



(B27) 



for a bond (ij) \\ ab. 

The respective system of equations of motion has a 
similar structure to that of Eq. (|B8|1 is: 
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( A\ 



-Bk 




*e k 

- n G (k) 
o 

5 k 



o 

-A k -n c (k) 



\ / ((b k \bt)) u \ 
-5* ((&k+q 

-A k -n c (k)J y ((6t_ k+Q | 6 t ))w J 



2,t | 



(B28) 



The Green functions can be now found from Eqs. (|B28p . 
They contain complete information about the bosonic 
correlation functions which appear in Eqs. (|B19p and 
(|B20p . They are obtained from the temperature Green's 
functions using the fluctuation-dissipation theorem^ 

k J 
k 

(B30) 

k J 



The values of the orbital correlation functions Eqs. (|B19|) 
and (|B20|) in the dimerized structure were used together 
with the respective spin correlation functions Eqs. (|B3[) 
and (|B3|) to obtain the self-consistent solution of Fig. [T3j 
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